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Abstract 

We study the effects of the perpendicular magnetic and Aharonov-Bohm (AB) flux fields on the 
energy levels of a two-dimensional (2D) Klein-Gordon (KG) particle subjects to equal scalar and 
vector pseudo-harmonic oscillator (PHO). We calculate the exact energy eigenvalues and normalized 
wave functions in terms of chemical potential parameter, magnetic field strength, AB flux field and 
magnetic quantum number by means of the Nikiforov-Uvarov (NU) method. The non-relativistic 
limit, PHO and harmonic oscillator solutions in the existence and absence of external fields are 
also obtained. 
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I. INTRODUCTION 



The pseudo-harmonic oscillator (PHO) potential is extensively used to describe the bound 
state of the interaction systems and has been applied for both classical and modern physics. 
It plays a basic role in chemical and molecular physics since it can be used to calculate the 
molecular vibration-rotation energy spectrum of linear and non-linear systems [1-3]. This 
potential is considered as an intermediate between harmonic oscillator (HO) and Morse-type 
potentials which are more realistic anharmonic potentials. In the non-relativistic quantum 
mechanics, the PHO is one of the exactly solvable potentials in the frame of the Schrodinger 
equation and has also been studied in one-dimensional (ID), two-dimensional (2D), three- 
dimensional (3D) and even in D-dimensional space. 

Therefore, this problem has attracted a great deal of interests in solving the Schrodinger 
and Klein-Gordon (KG) equations with the PHO interaction. The discrete (bound) and 
continuous (scattering) energy spectra of the PHO have been investigated by the SU(1,1) 
spectrum generating algebra [4]. The PHO potential has been obtained depending on the di- 
mension and the angular momentum and the zero-point energy in two-dimension was found 
to be minimum. The exact polynomial solution of the Schrodinger equation for PHO has 
been obtained in 3D space [5]. Recently, solutions of the Schrodinger and KG equations with 
PHO potential have been investigated [6-10]. A realization of the creation and annihilation 
(ladder) operators for the solution to the Schrodinger equation with a PHO in 2D was stud- 
ied in Refs. [11,12]. The operators satisfy the commutation relations of an SU(1,1) group. 
The exact solution of the Schrodinger equation with a PHO in an arbitrary dimension D 
was presented in Ref. [13]. The exact analytical solutions of the Schrodinger equation of 
the D-dimansional space for the PHO potential have been presented by means of the ansatz 
method and the energy eigenvalues were calculated from the eigenfunction ansatz [14]. The 
exact bound-state solutions of the KG and the Dirac equations with equal scalar and vector 
PHO potential have been obtained using the supersymmetric quantum mechanics, shape 
invariance and other alternative methods [15]. The bound-state solutions of the Dirac equa- 
tion with PHO have been obtained in the presence of spin and pseudospin symmetries [16]. 
In relativistic quantum mechanics, the treatment of this problem is relatively inadequate. 

Recently, several works have been achieved in this direction for which the solutions of non- 
relativistic and relativistic equations with different potentials for their bound and continuum 
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states were investigated. The approximate analytical solutions of the scattering states were 
derived for D-dimensional Schrodinger equation with modified Poschl- Teller, Eckart, modi- 
fied Morse and generalized Hulthen potentials were found for any angular momentum states 
using new approximation scheme to centrifugal term [17-20]. It was shown that the energy 
levels of the continuum states of modified Morse potential reduce to those of bound states at 
the poles of the scattering amplitude [19]. Using the supersymmetry and shape- invariance, 
exact bound state solutions of Schrodinger equation with Poschl- Teller double-ring-shaped 
Coulomb potential were presented [21]. Exact solution of the one-dimensional KG equation 
with mixed scalar and vector linear potentials were studied in the context of deformed quan- 
tum mechanics characterized by a finite minimal uncertainty in position using the momentum 
space representation [22]. New ring-shaped harmonic oscillator for spin-1/2 particles was 
studied and corresponding eigenfunctions and eigenvalues were obtained by using the Dirac 
equation with equal mixture of scalar and vector potentials [23]. Further, the approximate 
analytical solution of the Dirac equation with the Poschl- Teller potential has been presented 
for arbitrary spin-orbit quantum number in view of spin-symmetry [24]. 

Recently, the spectral properties in a 2D charged particle (electron or hole) confined by 
a PHO potential in the presence of external strong uniform magnetic field if along the 
z direction and Aharonov-Bohm (AB) flux field created by a solenoid have been studied. 
The Schrodinger equation is solved exactly for its bound states (energy spectrum and wave 
functions) [25,26]. So, it is natural that the relativistic effects for a charged particle under 
the action of this potential could become important, especially for a strong coupling. 

The aim of the present work is to study the exact analytical bound state energy eigen- 
values and normalized wave functions of the spinless relativistic equation with equal scalar 
and vector PHO interaction under the effect of external uniform magnetic field and AB flux 
field in the framework of the NU method [27,28]. The non- relativistic limit of our solution 
is obtained by making an appropriate mapping of parameters. Further, the KG-PHO and 
KG-HO special cases are also treated. 

The structure of this paper is as follows. We study the effect of external uniform magnetic 
and AB flux fields on a relativistic spinless particle (anti-particle) under equal scalar and 
vector PHO interaction in Section 2. We discuss some special cases in Section 3. Finally, 
we give our concluding remarks in Section 4. 
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II. KG EQUATION WITH PHO INTERACTION UNDER EXTERNAL FIELDS 



The KG equation is a wave equation mostly used in describing particle dynamics in 
relativistic quantum mechanics. Nonetheless, physically this equation describes a scalar 
particle (spin 0). Moreover, this wave equation, for free particles, is constructed using two 
objects: the four-vector linear momentum operator P M = ihd^ and the scalar rest mass 
M, allows one to introduce naturally two types of potential couplings. One is the gauge- 
invariant coupling to the four- vector potential {A^ (~^)} 3 which is introduced via the 
minimal substitution P^ — > P^ — gA^, where g is a real coupling parameter. The other, 
is an additional coupling to the space-time scalar potential S' con f(^) which is introduced 
by the substitution M — > M + S^onfC^)- The term "four- vector" and "scalar" refers to 
the corresponding unitary irreducible representation of the Poincare space-time symmetry 
group (the group of rotations and translations in (3 + l)-dimensional Minkowski space-time). 
Gauge invariance of the vector coupling allows for the freedom to fix the gauge (eliminating 
the non physical gauge modes) without altering the physical content of the problem. Many 
choose to simplify the solution of the problem by taking the space component of the vector 
potential to vanish (i.e., One may write the time-component of the four- vector potential 
as gA = V con f(r), then it ends up with two independent potential functions in the KG 
equation. These are the "vector" V COI1 {(~?) and the "scalar" S con f(~r^) potentials [29,30]. 

The free KG equation is written as 



(^ + MVfe G (^) = o. (i) 



(f + -2) 2 -{E- V coni (f)) 2 + (Mc 2 + S coni (f))' 



^,0) = O, (2) 



Moreover, the vector and scalar couplings mentioned above introduce potential interactions 
by mapping the free KG equation for a 2D single charged electron as 

where we used the transformation ~$ — > + -~A*. Further, the 2D cylindrical wave function 
■0(^,0) is defined as 

V>(f, 0) = -Le^(r), m = 0, ±1, ±2, . . . , (3) 
v 2ir 



where m is the magnetic quantum number. This type of coupling attracted a lot of attention 
in the literature due to the resulting simplification in the solution of the relativistic problem. 



The scalar-like potential coupling is added to the scalar mass so that in case when S con f(r) = 
iKonfC^), the KG equation could always be reduced to a Schrodinger-type second order 
differential equation as follows 
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(f + +2(E± Mc 2 ) Konf(^) + M 2 c 4 - E 2 



^(^,0) = O. (4) 



The potential Konf( ' ) is taken as the repulsive PHO potential [1-3,11]: 

V coni (rl = V ^-- r fJ, (5) 

where r and Vo are the zero point (effective radius) and the chemical potential, respectively. 
Further, the vector potential A may be represented as a sum of two terms, 1 = % + % 
such that ^ x 3t\ = if and ^ x ~% 2 = 0, where if = Bz is the applied magnetic field, 
and ^ 2 describes the additional Aharonov-Bohm (AB) flux field & AB created by a solenoid 
in cylindrical coordinates [31]. Hence, the vector potentials have the following azimuthal 
components [32] 

2 r ' 2vrr \ 2 2vrr J V y ! 

Hence, the bound-state solutions of the two cases in Eq. (4) are to be treated separately as 

follows. 



A. The bound states for positive energy 



The positive energy states (corresponding to S' con f( r ) = +V con f( r ) in the non-relativistic 
limit (taking E — Mc 2 — > E and (E + Mc 2 ) — > 2fxc 2 , where M is an effective mass and 
\E\ <C Mc 2 ) are solutions of 



^(^,0)=O, 



(7) 



where ip(~r^, 4>) stands for either ifj( + \~? , 0) or ip^ KG \~f^ , 0). This is the Schrodinger equation 
for the potential 2V con {(~?). Thus, the choice S con f(r) = +V con f(r) produces a nontrivial 
non-relativistic limit with a potential function 2V^ on f(~r^), and not V con f(~?). Accordingly, it 
would be natural to scale the potential term in Eq. (4) and Eq. (7) so that in the non- 
relativistic limit the interaction potential becomes V con {, not 2V con f. Thus, we need to recast 
Eq. (4) for the S(f) = V(f) as [29,33] 

c 2 (-ifr^ + -2) 2 + 2 (E + Mc 2 ) Vw(f) rl>(!>, 0) = (E 2 - M 2 c 4 ) 0), (8) 
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<K^,0) = O. (9) 



and in order to simplify Eq. (8) we introduce new parameters Ai = E + Mc 2 and A2 = 
E — Mc 2 so that it can be reduced to the form 

, 2 

Now, inserting Eqs. (3), (5) and (6) into the KG equation (9) and further introducing 
a change of variable s = r 2 , that maps r G (0, 00) to sG (0, 00), we obtain second-order 
differential equation satisfying the radial wave function g(s), 



(-iftf + ~ ( 1) 2 ~ \l (A 2 - Konf(^)) 



g"(s) + V(s) + i (-^s 2 + - L ) g ( s ) = 0, g(0) = and 0(00) = 0, (10) 



where we have employed the following abbreviation symbols 

1 



v 2 



h 2 c 2 



[Ai (A 2 + 2V ) - Mc 2 co c Hm'] , (11a) 



(3 2 = m' 2 + -^r*Vo\ u w! = m + rri = 1, 2, . . . , (lib) 
2 /M Wc A 2 1 VoAx 

7 = (^J (llc) 

with w! is a new quantum number. Here ^ = §ab/&o is taken as integer with the flux 
quantum $ — hc/e, u c = eB/Mc is the cyclotron frequency. Now, we use the basic ideas of 
the NU method [27] and the parametric NU derived in Ref. [28] to obtain energy spectrum 
equation 

v 2 = 2 (2n + 1 + 0) 7, n = 0, 1, 2, . . . . (12) 

where the constant parameters used in our calculations are displayed in Table 1. Inserting 
Eqs. (lla)-(llc), we finally arrive at the following transcendental energy formula, 



2 | 1 + 2n + l/m » + W 1 feV + w. 



ft 2 c 2 \ \ 2h h 2 



2^2 



"err. 







= ^ [Ai (A 2 + 2F ) - Mc 2 ^ c m'] , m' = 1, 2, . . . , (13) 

We may find solution to the above transcendental equation as E = E^q. In the non- 
relativistic limit when Ai — > 2Mc 2 and A 2 — > E, the above energy equation has a solution 
given by Eq. (39) of Ref . [25]. 

Using Eq. (38) of Ref. [28] and Table 1, we find the corresponding radial wave function 
g(r) as 

9(r) = C n , m r^e-^ 2 / 2 E (-n, \fi\ + l; 7 r 2 ) , (14) 



and hence the total 2D KG wave function (3) takes the explicit form 

*) = ^^^r*.^W>. (15) 

where (x) = ^^' F {—a, b + l;x) is the associated Laguerre polynomial and F(—a, b; x) 
is the confluent hypergeometric function. Note that the wave function (15) is finite and 
satisfying the standard asymptotic analysis for r = and r — > oo. 



B. The bound states for negative energy 



When S'confC' ) = —V con f(r), we need to follow same procedure of solution in subsection 
A and consider the solution given by (12) with the changes 

1 



2 . ~2 
V — > V 



h 2 & 



[A 2 (Ai + 2V Q ) - Mc 2 co c hm'] , 
1 



F^F = m' 2 + T^r z V X 2 , 



2 . ~2 

7 ^7 



~2h 



1 VqA 2 

/l 2 C 2 Tq 



(16a) 
(16b) 
(16c) 



Hence, the negative energy solution for antiparticle can be readily found as 



2 \2n + 1 + A/m' 2 + 



h 2 6 



■r 2 V \2 



~2h 



+ 



1 VbAg 

ft 2 C 2 Tq 



1 



fo 2 c 2 

and the 2D wave function is 



[A 2 (Ai + 2V ) - Mc 2 hoj c m'] , m! = 1, 2, . . . , 



(17) 



7r n + 







2tt 



(18) 



It should be noted that the negative energy states are free fields since under these conditions 
Eq. (4) can be rewritten as 



-hft + ffi +E \ w m (^,0)=o, 



2/i V" ' C "J ' ~J rn ' mv ' 1T ' ^ 
which is a simple free-interaction mode. Further, the set of parameters given in Eqs. (16a)- 
(16c) becomes 

_ [2ME M Wc ; ~ , . Mw c 

v — \ \ — — : — m', o — m , 7 = _ : ■ (20) 



ft 2 



h 



2h 
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Thus, Eq. (17) gives the following energy formula 



E ( J=(n + m> + ±)hw c , (21) 



and hence the wave function reads 



A • /e^ r2 L['"'» ( ^rA e im t (22) 
ix(n + m')\ n \ 2h J v ' 



III. DISCUSSIONS 



In this section we briefly study some special cases and relationship between our results 
and some other authors': 



A. Schrodinger-PHO problem under the effect of magnetic and AB flux fields 

In the non-relativistic limit (Ai — > 2M and A 2 — > E), Eq. (13) can be easily reduced into 
the form: 

E nm (£, (3) = hn(n+ -^—j + ^ c m' - 2V , (23a) 

Q = ^u 2 c + 4u 2 D , u D = ^2V /Mrl \m\ = Vm' 2 + a 2 , (23b) 

where a = k F r with k F = ^j2MVo/h 2 is the fermi wave vector of the electron [25]. Further, 
the wave function becomes 



with the given parameters 



h 2 \j\ 2h ) ' h 2 r 2 ' 

B. KG-PHO problem 

The energy spectrum of relativistic spinless particle in the absence of magnetic and AB 
flux fields has the form 



2U r V Q [l + 2n+\lm 2 + I = (A 2 + 2V ) y/X 1 , (25) 



r \ V h 2 



and is reduced to its non-relativistic limit: 



E = -2V 0+ ,l + 2„ + ^ + ^)^, (26) 

which is completely identical to Eq. (27) of Ref. [7]. The wave function can be expressed as 



with 



B = \jm^ + r ^\ l and C = ^r\/^V 



C. KG-HO problem 

The energy equation of the relativistic spinless particle subject to the HO field is 

n'hc\ r Zk - v^As = 0, ri = 1, 2, . . . , (28) 

where n' — 1 + |m| + 2n, n = 0, 1, 2, . . . , which is completely identical to Eq. (11) and Eq. 
(26) of Ref. [34] when one uses the notation k = 2V /r% = Mw%. Following Ref. [35], the 
Eq. (28) has three solutions, the only real solution giving energy is 

E nrn = - (Mc 2 + M 2 c 4 T~ 1/3 + T 1/3 ) , (29) 
o 

with 



T = 27kn' 2 h 2 c 2 - 8M 3 c 6 + 3n' 'Hey/ '3k (27 'kn' 2 h 2 c 2 - 16M 3 c 6 ). (30) 
The wave function takes the form 



= » xfn+ H)l rWg " l *' / ' £ ^ >(1 ^^' ° = <31) 



If one expands Eq. (28) as a series of A 2 , it becomes 

1 -o I 



/M 

n n = a/ — 



A2 + 4^ A ^32MV A2 + 0(A42; 



(32) 



and taking the first order of A2 by neglecting the higher order relativistic corrections, we 
finally arrive at the non-relativistic solution: 

E' nm = E nm - Mc 2 = (1 + 2n + \m\) = (1 + \m\ + 2n) hJ D , n = 0, 1, 2, . . . , (33) 

and wave function resembles the one given in (31) with D = -^\ / 2MV . 
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IV. CONCLUDING REMARKS 



In this work, we have obtained bound state energies and wave functions of the relativistic 
spinless particle subject to a PHO interaction and expressed in terms of external uniform 
magnetic field and AB flux field. We explored the solution of both positive (particle) and 
negative (anti-particle) KG energy states. The Schrodinger bound state solution is found as 
non-relativistic limit of the present model. It is noticed that the solution with equal mixture 
of scalar and vector potentials can be easily reduced into the well-known Schrodinger solution 
for a particle with an interaction potential field and a free field, respectively. We have also 
studied the bound-state solutions for some special cases including the non-relativistic limits 
(Schrodinger equation for PHO and HO under external magnetic and AB flux fields) and 
the KG equation for HO and PHO interactions. The results show that the splitting is not 
constant and dependent on the strength of the external magnetic field and AB flux field. 

Acknowledgments 

We thank the referee (s) for the invaluable suggestions and critics that have greatly helped 
in improving this paper. S.M. Ikhdair acknowledges the support of the Scientific and Tech- 
nological Research Council of Turkey. 

Appendix A: A Short Solution to KG Bound States 

Here, we seek to find a short solution to Eq. (9). By substituting Eqs. (3), (5) and 
(6) into Eq. (9), we obtain a 2D Schrodinger- type equation for the anharmonic oscillator 
V(r) = 7 2 r 2 satisfying the radial wave function R(r), 



where g(r) = r l l 2 R{r). For brevity, the solution of the above equation can be found through 
Eq. (17) of Ref. [5] by simply making the changes: 2nE' nL /h 2 — > v 2 , 2jjB 2 /h 2 — >■ 7 2 and 
2L + 1 -> 2/3 in Eq. (19) of Ref. [5] and then finally obtain our Eqs. (12) and (13). A first 
inspection on the asymptotic behavior of Eq. (Al), we find out that when r approaches 0, 
the radial wave function R(r) ~ r q , q = (3 + 1/2 > and when r — > oo, R(r) ~ exp (— r yr 2 /2) 




(Al) 
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and both solutions are satisfying the finiteness of the radial wave function g(r = 0) = and 
g(r — > oo) — > 0. Therefore, for the entire range r e (0, oo), we consider the general solution 
g(r) = r^exp (— r yr 2 /2) h(r), (3 > 0, where h(r) is the associated Laguerre polynomials, i.e., 
h(r) = L^n\br 2 ). With these behaviors, the 3D wave function (20) of Ref. [5] gives the 2D 
wave function as in Eq. (15) when one uses the relation: 




where N is the normalization constant. 
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TABLE I: Specific values of the constants in the solution of Eq. (20). 



Constants: (013 = case) 



£1 = 

s 1 


7 2 /4 




6 = 


z/2/4 

/ 


— 






"1 — 


1 

_L 


Q 2 = 


a 3 = 


04 = 05 = 


«6 = 


£1 = 7 2 /4 


Ctj = 




= -^ 2 /4 


a 8 = 


£3 = /? 2 /4 


Qg = 


«6 = 


7 2 /4 


"10 = 


= /3 + l 


an ~- 


= 7 




"12 = 


= /3/2 


"13 = 


--ll 


2 
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